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Abstract. When a physical system is put in contact with a very large thermal 
bath, it undergoes a dissipative (i.e., an apparently irreversible) process that leads 
to thermal equilibrium. This dynamical process can be described fully within 
quantum physics, involving only unitary, therefore reversible, maps. The infor- 
mation, initially present in the system, is not erased, but is diluted in the bath 
because of entanglement. Irreversibility may arise if, after quantum information 
has been thus diluted, some classical information is lost. This paper reviews a 
model for thermalization that displays these features. Two new analytical results 
are provided for the zero-temperature channels: a new quantitative measure of 
entanglement, and a study of irreversibility in the case where the lost classical 
information is the label of the particles in the bath. 

1 Introduction 

A well-established tenet of statistical physics says that an ensemble of physical systems in ther- 
mal equilibrium with a large reservoir ("canonical ensemble") will show a statistical behavior: 
a state of energy E will be occupied with a probability proportional to e~^^ where (3 = jt^, 
T is the temperature and fcs is Boltzmann's constant. Statistical physics, whose fruitfulness is 
beyond question, takes the existence of statistical ensembles as a starting point. The question 
about the origin, and even the meaning, of these statistics, is ultimately still unsettled. Re- 
cent remarkable developments have shown that ensembles are somehow generic in a kinematic 
sense: for instance, if one picks at random a state in the Hilbert space describing a large 
number N of particles, then, under suitable constraints, the state of much smaller sub-systems 
Pn = TrAT.nltZ/) (if"! shall almost always be close to a canonical state |1I2I3I4I5] . 

Here, we consider rather a dynamical process, the approach to thermal equilibrium, or 
thermalization, a typical example of dissipative phenomenorQ. It is defined as follows. A large 
number N of particles are already in the canonical state ^; one brings along a new particle, 
prepared in an arbitrary state p, possibly pure: thermalization is the process, at the end of 
which the new particle reaches arbitrarily close to the canonical state, while the bath is almost 



^ e-mail: physvOnus . edu . sg 

^ Some terminological issues to avoid confusion: I use dissipation as the phenomenological fact that 
energy, or information, which was initially concentrated in some physical system, flows through the 
evolution into the systems that interact with the first one. When a great number of degrees of freedom 
is involved, dissipation leads to apparent (and practical) irreversibility; but contemporary physics does 
not contemplate any fundamentally irreversible process. 
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unmodified. The naive description of tliis process, 

(1) 

is clearly non-unitary, because input orthogonal states are mapped on the same final state. 
For sure, non-unitary maps are allowed in quantum physics as descriptions of open systems: 
specifically, ([!]) is possible if the bath itself is coupled to an environment, to which it transfers all 
the information about p. Remarkably, quantum physics allows to find unitary maps that involve 
only the bath and the new particle, and that describe thermalization as well. This possibility 
is due to entanglement: the information, which was initially concentrated in p, is not lost in 
an external environment, but is encoded differently, being spread between the system and the 
bath in a coherent and reversible way. As a result, one may have 

p f>N — , ^ s^ch that TrATCT « ^ . (2) 

This means that the final state is not the exact thermal state of -I- 1 particles ^"^(^+1); 
however, when one picks any particle out, its partial state is as close as desired to the thermal 
state. In other words, Tr(A(j) will correspond to the expected thermal average for all single- 
particle observables A = jf^jp^ T^^=i -^k, with A]^ the operator acting as A on particle k and 
trivially on the others. In short: entanglement allows to construct a dissipative channel with 
fully reversible dynamics. The study of such channels is richer, and probably more satisfactory 
from the standpoint of physics, than the simple acceptance of the naive map (fT|). 

In this paper, I review the thermalizing channels that were presented a few years ago for 
a specific model of the system and the bath f6] . This approach to thermalization, that can be 
seen in the broader context of "quantum homogenization" 7 8, was proposed as a benchmark 
for exploring a quantitative link between entanglement and dissipation. I also present two new 
results for this model: a new measure of entanglement, and an analytical computation for a 
numerical result obtained in [7] for a model of irreversibility. 



2 The Model and the Results 
2.1 Definition of the Model 

The model is defined by the following assumptions on the kinematics and dynamics. 

Kinematics and Free Dynamics. The particle to be thermalized (the "system" S hereafter) 
is a two-level system, i.e. a qubit. The thermal bath is a reservoir composed of an arbitrary large 
number N of qubits, which interact with an external field but not with one another (a model 
featuring such interactions has been studied in Ref. [5]). The free Hamiltonian for the bath is 
therefore Hb — where h[k\ is the operator acting as ft. = —Eaz on the qubit k and 

trivially on the other qubits. We denote the projectors on the eigenstates of cr^ by Pq — |0)(0| 
and Pi = |1)(1|. The initial state of the system is arbitrary. The single-particle equilibrium 
state is 

e = e-'3VTr(e-'3'^) = pPo + qPi (3) 

with p = 5(1 + tanh(/?i?)) and q ^ 1 — P- We set i? > 0, so that |0) is the ground state, 
and p = 1 corresponds to T = 0. The initial state of the bath is the thermal state ps = 
g-/3ffB/Tr(e-^^«) = (C)®^. 

Interaction System-Bath. We consider a collision model, in which the system interacts se- 
quentially with the qubits in the bath one by one. Each step of this stroboscopic evolution is 
described by a hamiltonian H, or the corresponding evolution U ~ e'^, acting on (8) C^. 
Moreover, we consider that a qubit of the bath undergoes at most one interaction with the 
system (Fig. [T]), so the input state of the ancilla is always ^. The evolution of the system is 
finally described by the iteration of the completely-positive (CP) map defined as 



p(k+i) ^ xrs [U ip'^'^ (g> U^] = [p^'^^] 



(4) 
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Fig. 1. The coUisional model for thermalization: quantum circuit representation. The input (left) shows 
system qubit in an arbitrary state p, bath qubits in the thermal state ^. The unitary operations connect 
sequentially the system qubit with each of the qubits in the bath. At the output (right), the qubits 
are in a correlated state a. Thermalization is achieved if Trjvcr « ^ holds for all the single-qubit states, 
whatever the initial state p of the system qubit. 



2.2 Requirements for a Thermalizing Channel 

We want to find all the two-qubit unitary operations that define a thermalizing channel in our 
scenario. Four requirements are suggested by physics: 

First Requirement. U should not depend on temperature (i.e. onp), but may depend on the 
label z of the local free Hamiltonian. The requirement is motivated by the fact that z is a locally 
available information, while temperature is a characteristic of an ensemble. Also, we want to 
characterize interactions that thermalize at all temperatures, avoiding possible pathological 
examples which would define the correct physics only for a specific value of this parameter. So 
from now on we write Uz- 

Second Requirement. If the system is prepared in the local equilibrium state ^, nothing 
should happen, because the thermalization is already achieved. Formally: 

Uz{(,®£,)Ul^^®i for ^=pPo + qPi. (5) 

Note that this is stronger than requiring ^ to be a fixed point of T^, because we ask that the 
qubit of the bath is unchanged too. Because of the first requirement, ^ should hold for all p. 

Third Requirement. For any input state p of the system, the iteration of the map leads 
to thermalization: 

^T^lp] yp. (6) 

Fourth Requirement. For any input state p of the system, the fluctuations introduced in the 
state of the bath are small. 



2.3 Solution: All Thermalizing Channels for the Model 

The combination of the first and the second requirements imply that the subspaces Pq ^ Pq, 
Pi (8) Pi and Pq^Pi + Pi^ Pa must be invariant under the action of Uz . To prove this assertion, 
we notice that on the l.h.s. of ([5]) the term Uz Pq ^ Pq Ul appears with the weight , the term 
Uz {Pq (E)Pi + Pi(E) Pq) Ul with the weight p(l — p), and the term Uz Pi ® Pi Ul with the weight 
(1 — p)^. Since we want condition © to hold for all p, the three subspaces must be separately 
invariant. Thus the first and the second requirement restrict U to take the form 



|00) — > e'^«\00) 

jll) — > e^'^Mll 
|01) 
|10) 



cos<?!)|01) +e*'^=^sin0|lO)] '^''^ 
cos</)|10) -e-''^=^sin0|Ol)] . 
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Remarkably, each of these unitaries defines a thermalizing channel (in other words, the third 
requirement is automatically fulfilled when enforcing the first and the second). To prove this 
assertion, we compute explicitly p'"' = T^[p] as a function of the parameters of the initial state 
p. Let's write 

pin) ^d(")po + (l-d("))Pi + /c(")|0)(l| + fc(")*|l)(0|. (8) 

Inserting the explicit form ([7]) for Uz into (jl]) , we find that the effect of the map is given by 
^ (i(")cos2(/) + psin2^ and = cos0Afc(") with A = pe*(xo-X3) _^ ^gi(x2-xi) (note 

that |A| < 1). A straightforward iteration gives c?^"^ and fc^"^ as a function of the parameters 
d*^"' and fc*^"' of the initial state p: 

= [1 - (cos 0)2"] p + (cos 0)2" , (9) 
fc(«) = (A COS0)". (10) 

Thus, whenever <j) ^ 0, the iteration of yields d^"-' p and fc*^"^ — > 0, i.e. thermalization ([6]). 
Finally, the fourth requirement is fulfilled by taking small enough. In fact, it can be verified 
that the fidelity of the state cr of a bath qubit after the interaction with respect to the thermal 

state satisfies = Tr (^^ia^i)^ > cos0 for all temperatures and for all input states of the 

system. 

Of the six free parameters defining the most general Uz ©, only (j> and the differences 
Xo — X3 ^nd X2 — Xi define the CP-map T^. The reason is that there are two symmetries of the 
physical process of thermalization. The first one is the usual choice of a global phase; the second 
is the freedom of choosing the global phases of |0) and |1) for qubits in the bath, both before 
and after the interaction with the system. Mathematically, if Uz defines a CP-map T^, then, 
writing u{x) = Pq + e^^Pi a rotation along the z axis, e*-^ [l (8) u(a)] Uz [l (8) u(/3)] defines the 
same CP-map for all choices of a and /?. So much for the three parameters which play strictly 
no role. In addition, the phase of A is a measurable parameter but does not add any insight 
on the thermalization process. In fact, this parameter is associated to a rotation around the z 
axis, or in other words, to a redefinition of the x and y axes of the Bloch sphere in the plane 
perpendicular to z, for the system qubit. Since there is nothing, in the physics of the model, 
that would single out a specific axis in the (x, y) plane, we can study thermalization for any 
specific choice of the phase of A. 

All in all, the physics of thermalization in this model can be studied on a two-parameter 
family of unitary transformations. A convenient choice 

|00) — > e*^|00) 

|oi) _^cos0|Ol)+isin0|lO) ^ ^ 

1 10) — > cos 01 10) +isin0|Ol) 

which can be written as e^^/^e*^^"^^^' with 

6*) = ^ [0 [cTa, (g>a^+ay(g> ay) + doz ® cr^] • (12) 



2.4 The Physics of Thermalizing Channels 

In spite of the extreme simplicity of the model, the thermalizing channels display a rich physics. 
Here is a review of these features. For the known results, the demonstration is given in Refs 
|6I7| ; the new results will be demonstrated in Section [3] below. 



^ Here the sign of 6 is the opposite than in Ref. [B]. 



Will be inserted by the editor 



5 



The partial swap. The particular choice 9 = (p makes Vz{(t), (p) = V{(l)) independent of z: this 
is the only interaction, in this model, that is independent not only of temperature, but also on 
the parameter of the local Hamiltonian. This unitary can be written as 

V{<j)) — cos + i sin (pUswap (13) 

where Uswap — ^i'^/'^) is the operation that swaps the state of the two qubits: |V'i)|'02) ~* 
|V'2)|V'i)- The mechanism for thermalization is quite intuitive in this case. With a channel made 
of pure swaps, the system would be in state ^ after one single interaction, but this would 
introduce a large fluctuation in the bath. The partial swap, in the meaningful limit coscj) ~ 1, is 
the gentle version of this process: after many steps, the state of the system has been "swapped" 
into the bath, but the fluctuations in the bath itself are small. 

Dissipation and Decoherence. The dynamics of thermalization is described in ([9]) and (jlOp . 
with A = pe*^ + qe~^^ . The parameter 9 appears only in the dynamics of the off-diagonal term, 
hence is entirely related to decoherence. The re-equilibration of populations, i.e. dissipation, is 
governed by (f> alone0. Interestingly, dissipation and decoherence are decoupled in this model, 
in the sense that T4(0, d) = K(0, 0)T4(0, 9) = Vj(0, 6I)K((/), 0). Also, one can verify that 

V,{^,9) ^V,{0,9~(b)V{(b) : (14) 

all the thermalizing unitaries in this model can be seen as the application of the partial swap 
followed by an additional reduction of coherence. 

Continuous-time limit. In order to pass from the discrete dynamics indexed by n to a 
continuous- time dynamics with parameter i, one sets n — t/ro, and lets the interaction time 

To go to zero together with (p and 9, keeping constant the dissipation rate ^ = ^ and the 

phase fluctuation rate = ^^j. One finds that in the continuous-time limit, the processes of 

dissipation ^ and decoherence (fTU]) are exactly exponential: 

d{t) = e-*^^'d{0) + e~''^')p , (15) 
|fc|(0 = e-*/^^|fc|(0) (16) 
with ^ — +pq7^. For = or at zero temperature, the bound Ti > (see e.g. |10) . 

p. 120) is saturated. Note also that it is possible to relate this dynamics to a master equation 
of the Lindblad type [8]. 

Dissipation and Entanglement. As we mentioned in the introduction, this model was pro- 
posed as a benchmark for studying the link between entanglement and dissipation. This link 
has been explored along several directions. An especially strong link is provided by the following 
observation: we have seen that any thermalizing unitary is equivalent to one of the V{(j), 9) up 
to local unitaries (LU), and it can be shown that all the V{(t)^9) are inequivalent under LU 
|6|llj . Thus (j) and 9 (or any equivalent choice of two parameters) are necessary and sufficient 
to define all the properties of entanglement; but in turn, these parameters uniquely define the 
relaxation times Ti and T2. So, for the model under study, the relaxation times are directly 
related to entanglement. 

In Refs [617112] several kinds of entanglement have been computed: the largest amount of 
two-qubit entanglement that a single application of V{4), 9) can generate, the entanglement 
of a given qubit versus all the others etc. These are measures of bipartite entanglement, the 
only ones which were available at the moment of writing. Computable measures of multipartite 
entanglement have been proposed since [13 14115] . According to one of these measures, after 
interaction of the system with n qubits in the bath, an amount of entanglemen10 



1 _ ^2(n-|-l) - 

£:„+i = 2|ci| ^/ (1 - c2«) ( 1 - |ciP ) , (17) 



^ Note that 4> has to appear in both the diagonal and ofT-diagonal term, as is indeed the case: any 
re-equilibration of populations must be balanced by a sufficient amount of variation in the coherence, 
in order to avoid evolution into non- physical states (matrices with negative eigenvalues). 

* For comparison, the same measure of entanglement applied to the (n + l)-qubit GHZ state gives 
£r,+i{GHZ) = y2Vl-2-". 
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with c = COS0, is generated at T = for an arbitrary pure initial state co|0) +ci|l) (subsection 
I3.ip . This expression is monotonically increasing in n; for any given n, the maximal amount of 
entanglement is generated by the input state |1), which is indeed the farthest from equilibrium. 
For this case ci — 1, £oo — 2yj=f V2 in the limit of weak interactions. 

Irreversibility. The thermalization process that we described is certainly reversible, since it 
is described by the unitary operation U = Us,n---Us,2Us,i- To study irreversibility, one must 
add some lack of knowledge. In our context, the most natural way of doing this is to suppose 
that the labels of the qubits in the bath are randomly permuted before reversing the evolution 
U, and we don't control this permutation. Then the reversed evolution reads — UU which 
is equal to the identity if and only if 77 = 1. Note that 77 is unitary, and consequently so is Ut^: 
apparent irreversibility is due to the fact that we reverse only U, because we are not supposed 
to control 7J. In Ref. [7], we have provided numerical results for this effect. In subsection 13.21 
below, we compute analytically the average fidelity F of the reconstructed state with respect 
to input state of the system qubit, for the case T = and for the specific evolution F(0, 6 — 0). 
Taking the excited state |1) as input, one has 

F.i.'Mi^ (1.) 

with c = cos ((ii, for n large enough; the expression for an arbitrary input state is Eq. (jSSp . From 
these analytical results, one learns that the expected decrease of the average fidelity is much 
slower than the decrease in 1/n! of the probability of retrieving the initial state exactly. 



3 Demonstration of the New Results 

This is a technical section, devoted to the demonstration of the two new results. The common 
starting point is the form of the (n + l)-qubit state obtained after interaction of the system 
qubit with n qubits in the bath. We consider an arbitrary pure input state 1-0) = co|0) + Ci|l); 
at T = 0, the equilibrium state is ^ = Pq pure as well, so the initial state is \^in) — \'4') g\Q ■ . .0) g 
and the state stays pure under the evolution. After the system has interacted with n qubits in 
the bath, the state reads 



n 

Z^nltf'm) =coe^"'|0)s|0")B + ci[c"|l)s|0")B+«s|0)s(Ec'"'e'^""'^'|lfe> 



k=l 



(19) 



where |0")^ and llfe)^ are n-qubit product states of all |0), respectively of |1) for qubit k and 
|0) for the others. Through all this section, we use c = cos(j), s = sui(j) and B — {1, ...,n}. 



3.1 Amount of multipartite entanglement 

We want to compute the multipartite entanglement in the [n + l)-qubit pure state (|19p . The 
measure of multipartite entanglement that we are going to compute is Eq. (6) of 13J, Eq. (88) 
of [H]: 

Sn+i = 2^-^ ^2"+i - 2 - Sn+l (20) 

where is Sn+i is the sum of Tr(p|) over all partial traces. Since the state is pure, the spectral 
properties of ps,b and of PB\b are identical for all subset h of B\ therefore we must compute 
Sn+i = 2 X^b^iiB TrpI J. Now, from one has that all the ps,b are rank two states of the 
form \^s,b){'^s,b\ + a;f,|0...0)(0...0| with \'fs.b) a non-normalized state whose precise form is not 
important here, and with 



cr, = l-||^S,.||' = |ci|V^c2('=-i) 

k^b 



(21) 
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So the quantity to be computed is 

Sn+i^2Y,[{l~xtf+xl] = 2"+i-2-4[^Xfc-^x^]. (22) 

b^B b^B b^B 

To compute the first sum, notice that f^^pp- J2b^B = J2b^B J2k<^b c^^''^^^ = Y^k=i c^'-'^^^^^fc 
where is the number of sets b such that k ^ b. Clearly ~ 2"^^ for all k; therefore we are 
left with a geometrical series and finally 

= 2"-i|cini-c2") . (23) 

b^B 



The second sum that we have to compute for ((22|) is 

b^B b^B k^b b^B k,k' ^b\k^k' 

1 1 _ 4n 

Proceeding exactly as above, one obtains Ii = 2"^ ^y^-S^^. Similarly, 

, f-Y- 2ik+k'-2)^ ,„_i e^(l-c^("-i))(l-c^") 
fc=i fc'^fc ^ ' 

because Nkk' = 2"^^ for all k and fc' 7^ fc, and we have applied to x = c? the generic formula0 



^ 2x(l-x"-i)(l-x") 
(1 - a;)2(l + 

Inserting everything back into (|^. then into one obtains pT)) 



t^.^. (l-.)^(l+x) 



3.2 Irreversibility via random permutations in the reservoir 

We start again from state (|19p. After the interaction with n qubits of the bath, we keep the 
system qubit and apply a random permutation tt : B — > B to the indices of the qubits in the 
bath. At this point, the state reads 

n 

™„|tf'„) =coe^"^|0)5|0")5 + ci[c"|l)5|0")5+z.|0)^(5]c*W-ie^("-*W)^|U)^)](27) 

fc=i 

with TT = TT^^. On this state we apply the reverse evolution, at the end of which we find 

UlnU,,\^,r.) = (co|0) +ciMl))^|0")5 + ciVl-|/.P|0)^yafc|lfc)5 (28) 

k 

where we don't need the explicit form of the ak, but we need 

n 
k=l 

Let pCn-.'r:"^) = Trs 0...0) (V-, 0...0|Z^i] be the reconstructed state of the system qubit: 

its fidelity with respect to the initial state is 

FiTT) EE = |coP + |cip [I/.P + 2|cop(Re/. - l/.H] . (30) 



^ This formula can be re-derived from the usual sum of geometric series by noticing that 
E.^^. equal to 2 E::: E:-=.+i . or ahematively to ( ^'f " 
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To quantify irreversibility, we must now average over all possible permutations. At this stage, 
we simplify the problem by focusing only on the unitary operation V{(j),6 = 0), so that /^r 
becomes rea0. To find F — Fij^) we have to compute 

1 " 

= c^" + s^i- y y c'=+*(fc)-2 ^ + s^h , (31) 

TT A:=l 

1 " 2 

Jl = + 2c2",s2/i + 5] ( ^ c'^-+*(^)-2^ = C^" + 2c2"s2/^ + ^4^^ (32) 

TT fc=l 

The pattern for these calculations is quite similar to the one we used in subsection |3Tj Consider 
first /i: observing that there are (n — 1)! permutations such that 7r(fc) — j, one can decouple 
the two sums and write 



k=i ]=i ^ ' 

Similarly, using in particular (|26p . we obtain 

^ ^ ^ ^2(fe+S-(fe)-2) + _ ^ ^ ^fe+fc'+#(fc)+#(fc')-4 

~ n 1 - c2 j n(n - 1) 1^ (l-c)2(l + c) 
So finally the average fidelity for the reconstructed state is 



(34) 



^=|coP + |ciP [/2 + 2|coP(A-/2)] « |coP + |cip [s% + 2|coP(s2/i-s4/2)] . (35) 

Inserting the expressions above, one finds that the convergence F — > |coP is indeed present as 
expected, but is very slow (remember that c = cos0 is close to 1 in the meaningful limit). The 
expression (ITS)) given above is obtained by setting cq = and c" w 0. 



4 Conclusions and Open Questions 

In conclusion, we have reviewed a coUisional model for the dynamical process of thermalization 
(or more generally, homogenization) . In spite of its simplicity, the model exhibits a rich variety 
of interesting physical features, e.g. decoupling of dissipation and decoherence and exponential 
decays in time for both (Ti and T-i). Moreover, this kind of model provides a benchmark to study 
quantitative links between the parameters of thermalization and the entanglement generated by 
the evolution. This relation has been studied in this paper using a new measure of multipartite 
entanglement. Also as expected, if classical information (here, the labels of the particles in the 
bath) is lost after the entanglement has been distributed, the channel becomes irreversible. 
Here, this irreversibility has been quantified analytically for an example of channel. 
Here is a list of open questions related to thermalization: 

— Staying within the model, the two new results presented here are partial, and suggest them- 
selves lines for further study. First of all, both results have been derived only for the case 
T = 0, because in this case the state of all the qubits is pure. In addition, the measure of 
multiparticle entanglement ()20p has no clear-cut operational meaning and, as it turned out, 

^ This is not the only simplification: if 6* 7^ 0, expression ([29} contains both 7f(fc — 1) and 7f(/c). 

^ It is known that this and similar measures of multipartite entanglement may detect entanglement 
beyond the bipartite case [15]. We note also that, contrary to what is stated in [^, the result obtained 
in [12] does not exclude the existence of multipartite entanglement in our thermalizing channels. 
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does not grasp all the physics of entanglement in the model, because the resulting quantity 
(|17p is independent of 9. Similarly, the irreversibility has been studied only for the channel 
with — Q, and in the case where only the bath qubits are permuted but the system qubit 
is known. 

— The studied model has some remarkable features: e.g., the fact that the third requirement 
is automatically implied by the two others; the fact that (up to local unitaries) thermalizing 
channels are characterized by two parameters in one-to-one correspondence with the relax- 
ation times Ti and T2; etc. Are these generic features of coUisional thermalizing channels, 
or just artefacts of the choice of deahng with qubits? 

— We discussed in the text the intuitive character of the partial swap. This operation is only 
possible because the system is of the same dimensionality as the particles in the bath. What 
happens if one wants to thermalize a d-dimensional quantum system with a bath of (say) 
qubits? Can the qubits be used one by one, or must one take them in bunches? 

I am very grateful to the organizers of the 382. Wilhelm und Else Heraeus Seminar (Bad Honnef, 
8-10.01.2007) for the invitation to attend this very stimulating workshop and give a talk there. This 
article benefits from discussions with Andreas Buchleitner, Daniel Burgarth, Oscar Dahlsten, Dai Li 
and Peter ZoUer. 
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